Directed Ising type dynamic preroughening transition in one dimensional interfaces 
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We present a realization of directed Ising (DI) type dynamic absorbing state phase transitions in 
the context of one- dimensional interfaces, such as the relaxation of a step on a vicinal surface. Under 
the restriction that particle deposition and evaporation can only take place near existing kinks, the 
interface relaxes into one of three steady states: rough, perfectly ordered flat (OF) without kinks, 
or disordered flat (DOF) with randomly placed kinks but in perfect up-down alternating order. A 
DI type dynamic preroughening transition takes place between the OF and DOF phases. At this 
critical point the asymptotic time evolution is controlled not only by the DI exponents but also by 
the initial condition. Information about the correlations in the initial state persists and changes the 
critical exponents. 

PACS numbers: 64.60.-i, 02.50.-r, 05.70.Ln, 82.65. Jv 



I. INTRODUCTION 

Absorbing type dynamic phase transitions are the fo- 
cus of extensive research [p]-|l2[] . These transitions occur 
in dynamic processes with trapped states. At the ab- 
sorbing side of the phase transition, the system evolves 
into one specific microscopic state, a so-called absorbing 
state, out of which it can not escape. At the active side of 
the phase transition, the system manages to avoid such 
traps. An ensemble has a finite probability to stay alive. 

Two distinct types of absorbing phase transitions have 
been identified in one dimension: the directed perco- 
lation (DP) and directed Ising (DI) universality class. 
DP type dynamic critical behavior has been found in, 
e.g., Sehlogl's first model for contact processes pair 
contact processes [B, and branching annihilating ran- 
dom walks (BAW) with an odd- number of offspring ||. 
DI type dynamic critical behavior has been found in, 
e.g., probabilistic cellular automata ||, nonequilibrium 
kinetic Ising type models |t| , interacting monomer-dimer 
models [p|, and BAW models with an even-number of off- 
spring MM ■ These models describe a wide range of phe- 
nomena, in particular epidemic spreading and catalytic 
chemical reactions. 

In analogy with equilibrium phase transitions, it is be- 
lieved that dynamic universality classes are determined 
by the degeneracy and symmetries of the absorbing states 
and by the symmetry p ro p erties of interfaces (domain 
walls) between them |Hj-|l2|, DP type transitions involve 
typically only a single absorbing state or a set of absorb- 
ing states with one of them dynamically more prominent 
in a coarse grained sense than the others |Q,10|,^2|. DI 
type transitions involve two equivalent adsorbing states 
or two equivalent classes of absorbing states [ l|| . 

Counting the number of absorbing states is not suf- 



ficient. The degeneracy of the absorbing states can be 
obscured by the formulation of the dynamic rule. For 
example, in BAW dynamics involving only one species of 
particles, the "empty" absorbing state seems not to be 
degenerate, but the transition belongs to the DI univer- 
sality class if the dynamics conserves the particle number 
modulo 2 ||[| . In that case particles can be reinterpreted 
as domain walls and we can color the domains at oppo- 
site sides of the walls alternatingly with two colors. One 
of the two colors dies out in the adsorbing state. This is 
reminiscent of the domain wall formulation of the equilib- 
rium Ising model, where the existence of two coexisting 
phases is also obscured. 

One appealing scheme of classifying absorbing phase 
transitions is by association with conventional equilib- 
rium phase transitions. For example, consider the one- 
dimensional Ising model with single spin- flip dynamics. 
In terms of domain walls this involves the following pro- 
cesses: spontaneous creation of domain wall pairs — » 2 A 
with probability p sp (the Boltzmann weight) , hopping of 
single domain walls with probability ph = 1, and pair 
annihilation 2 A — > with probability p a = 1. The 
stationary state of this model is the equilibrium state, 
which is obviously "disordered" (active) in one dimen- 
sion. However, the model evolves always into one of the 
two perfectly ordered (absorbing) states if we disallow the 
spontaneous creation of domain walls, p sp — 0, because 
then the domain wall density can only decrease. The 
DI transition comes into play when branching processes 
are allowed; i.e., the creation of domain wall pairs in the 
vicinity of existing domain walls, A — > 3A. Branching 
can keep the system alive, while the two perfectly or- 
dered Ising ground states still remain absorbing states. 
The resulting model is known as the BAW model with 
two offspring. 
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The same line of reasoning associates a distinct ab- 
sorbing type dynamic universality class with each con- 
ventional equilibrium universality class, like the familiar 
g-state Potts and g-state clock universality classes. How- 
ever, none of the models studied thus far with higher 
symmetries than the Ising model has an absorbing phase 
transition. Dynamic processes with more than 2 equiv- 
alent absorbing states (q > 3) appear to be always ac- 
tive |Tl],[l4|]. Cardy and Tauber provide a possible ex- 
planation for this |p.5| . Additional systematic numerical 
studies are needed to settle this issue. 

In this paper we introduce an application of absorb- 
ing phase transitions to the dynamics of one-dimensional 
solid-gas interfaces. We find a one-dimensional dy- 
namic preroughcning (PR) transition from an ordered 
flat (absorbing) phase to a disordered flat (active) phase. 
Consider the one-dimensional restricted solid-on-solid 
(RSOS) model description of crystal surfaces. The height 
difference between nearest neighboring columns of parti- 
cles is or ±1, i.e., only steps (kinks) with single atomic 
heights are allowed. A conventional dynamic Monte 
Carlo type equilibration process where single particles 
adsorb on or desorb from the surface with equal proba- 
bility leads to a rough equilibrium surface with Edwards- 
Wilkinson (EW) type |l6| scaling properties. In analogy 
with the above directed Ising model discussion, we trans- 
form the perfectly ordered flat state into an absorbing 
state by disallowing adsorption and desorption of parti- 
cles at flat segments of the surface (spontaneous creation 
of step pairs). The only processes allowed are: hop- 
ping/pair annihilation of steps by adsorbing or desorb- 
ing a particle at step edges and branching of steps by 
adsorption or desorption at the next nearest neighbour 
sites near existing steps. The branching processes fall 
into two classes: the ones that preserve local flatness and 
those that create local roughness. We give them inde- 
pendent probabilities. 

This model describes the evolution of the one dimen- 
sional interface from any initial random rough configu- 
ration to three types of stationary states: a rough phase 
with dynamic exponent z — 2 (EW type), a perfectly 
ordered flat (OF) absorbing phase without steps, and 
an active disordered flat (DOF) phase. DOF phases are 
well known in equilibrium surfaces. They represent step 
liquids where the surface remains flat due to long-range 
step-up step-down alternating order but steps are placed 
randomly [17|,[l8). Our dynamic DOF stationary state 
has no defects, i.e., the up-down alternating order is per- 
fect and the the steps are placed randomly. 

Although the dynamic DOF phase is active, it has an 
absorbing state type property. The OF absorbing state 
and the DOF active phase at the other side of the tran- 
sition have in common that both "ground states" lack 
thermodynamic defects (both are at their fixed point in 
the sense of renormalization theory). 

Equilibrium preroughcning transitions belong to the 
Ashkin- Teller universality class 1 17 . They involve two 
coupled Ising type order parameters. They are non-zero 



at opposite sides of the transition. A perfectly DOF ini- 
tial configuration has a half integer average surface height 
n + ~ and evolves at the OF side of the PR transition into 
a perfectly flat ordered absorbing state with n or n + 1 
integer surface height. This OF- type two- fold degener- 
acy is the conventional spontaneous symmetry breaking 
associated with DI type transitions. However, an Ising 
type degeneracy exists also at the DOF side of the PR 
transition. An initial state with an integer average sur- 
face height n, evolves into a stationary DOF state with 
surface height n + i or n — |. We call this the DOF- type 
degeneracy. In analogy with equilibrium transitions one 
might expect therefore that the dynamic PR transition 
be described by 2 coupled DI transitions. This turns out 
not to be the case. 

The dynamic PR transition belongs to the (single) DI 
universality class. The two types of spontaneous symme- 
try breaking are not on an equal footing. A perfect DOF 
initial state can decay into the OF state, but a perfect 
OF state is frozen forever (even at the DOF side of the 
PR transition). The DOF type degeneracy does not cre- 
ate additional DI critical fluctuations. But it still affects 
the scaling behaviour. Configurations with perfect DOF 
order form an invariant subspace (which includes the OF 
ordered state). As shown below, this leads to a peculiar 
initial configuration dependence of the asymptotic decay 
of the step densities (the kinks in the one-dimensional 
interface). 

In section 2, we introduce our model in detail. It is 
instructive to interpret it not only as a model for surface 
relaxation, but also for surface catalysis. In its latter 
reincarnation, the process is a two-species generalization 
of the BAW model with two offspring. The up and down 
steps (kinks) represent two species of particles, A and B. 

In section 3 we limit ourselves to configurations with 
perfect ABAB alternating order. These DOF type con- 
figurations form a dynamical invariant subspace. The 
numerical results presented in section 3 confirm that our 
PR transition belongs to the DI universality class. 

In section 4, we discuss the crossover scaling proper- 
ties of the DI critical point into the rough phase. The 
rough phase has conventional EW type scaling proper- 
ties. Recall that we disallow spontaneous adsorption and 
desorption from flat surface segments; surface roughness 
can only be created and maintained by branching. Ap- 
parently, this restriction does not alter the scaling prop- 
erties of the rough phase. 

The DOF-type Ising degeneracy shows up in the evo- 
lution of arbitrary initial states. The kink densities scale 
in time in an unusual way, not only at the PR transition, 
but also everywhere in the OF and DOF phases. They 
depend strongly on the initial conditions; whether the 
initial configuration is flat or rough, and on the correla- 
tions in such initial rough states. Conventional wisdom 
tells us that the long time scaling of dynamic processes 
depends only on the dynamic exponent and the station- 
ary state exponents of a process. In our case, critical 
exponents vary with the initial conditions. We present 
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an analytical scaling theory for this in section 5 and nu- 
merical results in section 6. 



II. THE MODEL 

Consider a one-dimensional lattice. Each site is vacant 
or occupied by at most one A or one B type particle. 
In the surface catalysis interpretation, A and B repre- 
sent two types of particles. In the surface formulation, 
they represent up and down steps. Configurations evolve 
in time according to the following dynamic rule. First 
choose a site at random. If the site is empty, nothing 
happens. If the site is occupied, the particle can hop to 
its nearest-neighbor sites with probability p. If it lands 
on top of an existing particle of the opposite kind, the 
AB pair annihilates immediately. The move is rejected 
if the particle would land on top of a particle of the same 
type. 

Besides hopping, each particle can also branch into 
three particles by the creation of an AB pair. Branch- 
ing comes in two distinct flavours; the one that preserves 
local ABAB order and the one that breaks it. Order 
preserving branching (OPB) creates local DOF step-up 
step-down order in the surface. Order breaking branch- 
ing (OBB) creates local roughness (see Fig. [l]). These 
branching processes occur with probability q and r, re- 
spectively. Branching could lead into double occupancy 
of a site. AB pairs at the same site annihilate immedi- 
ately. The branching event is rejected if it would result 
in two particles of the same type at any site. We require 
that the chosen particle must attempt to hop or branch. 
This implies that p + q + r = 1. q = q/p and f = r/p are 
our choices for independent parameters. 

We need to distinguish also between so-called dynamic 
and static branching. The two new particles, A and B, 
created by each branching process may be placed in two 
different ways. The center of mass is stationary (static 
branching) or moves (dynamic branching) fTsf| . In the 
latter the new particles are placed both to the left or 
both to the right of the parent particle with equal prob- 
ability. In the interface representation of the model, dy- 
namic branching corresponds to single particle adsorp- 
tion/desorption. Static branching represents two-particle 
events (see Fig. |l|). So dynamic branching is more nat- 
ural. These details do not affect the universality of the 
phase transition. But they can change dramatically the 
location of the critical point jl9) . It is well known that in 
the BAW model with two offspring the stationary state 
is always absorbing. Static branching does not create 
enough activity to destabilize the absorbing phase poj . 
We applied both type of branching in our simulations. 
The universal scaling properties do not change. Here 
only our results for dynamic branching are presented. 

Each particle configuration of the surface catalysis 
model maps onto a surface height configuration of the 
RSOS type surface growth model. The A and B par- 



ticles represent up and down steps of unit height. The 
above dynamic rules for the A and B type particles trans- 
late into placement or removal of a single surface atom 
near existing steps (Fig. [I]). Adsorptions and desorptions 
have equal probability. The surface does not advance nor 
retreat on average. Flat segments of the surface are in- 
active, which makes the perfectly ordered flat state an 
absorbing state. 

The structure of the phase diagram is shown in Fig. |^. 
The numerical details will be presented in section 3 and 
4. Here we want to point out the general features. 

OPB type branching preserves the average surface 
height, while OBB type branching creates local rough- 
ness. This suggests two distinct types of active phases. In 
the region of the phase diagram where OPB type branch- 
ing is dominant, the steps prefer up-down alternating or- 
der. The surface likes to remain flat on average with ran- 
domly placed steps. This is a one-dimensional dynamical 
version of the DOF phase known from two-dimensional 
equilibrium crystal surfaces Jl7],[l8| ■ In the region of the 
phase diagram where OBB type branching is dominant, 
the active phase lacks step up-down alternating order and 
is therefore probably rough. 

In the presence of OBB type branching (f > 0), the 
stationary surface state turns out to be always rough 
(see Fig. |]). Along the f = line, the surface is flat; 
in the OF absorbing phase for q < q c and the DOF ac- 
tive phase for q > q c . Spontaneous order is difficult to 
maintain in one dimension. This explains why the DOF 
phase is limited to the f = subspace. It would have 
been nice, but a surprise if the OF phase extended into 
f > 0. Equilibrium roughening phase transitions can be 
viewed as q — » oo limits of the g-state clock model. Ab- 
sorbing phase transitions with q > 3 are unknown and 
may not exist at all as mentioned in section 1 |21|. 

The transition point q c in Fig. [2] represents a onc- 
dimensional dynamic analogue of equilibrium prerough- 
ening transitions. The disordered flat phase is active (a 
step liquid phase), but maintains perfect step-up step- 
down alternating order. A and B type particles move 
around like in a liquid, but remain perfectly ABAB or- 
dered. The distances between neighboring A and B par- 
ticles are randomly distributed. 

To determine the scaling properties of this transition 
(section 3) we measure the step (particle) density p(t) as 
a function of time t, which is the sum of the two single 
particle densities, pA{t) + pB(t)- The difference between 
them, pA(t)—pB(t), describes the global tilt of the surface 
and is preserved by the dynamic rule. The step density 
p in the stationary state vanishes in the absorbing phase 
and remains non-zero in the two active (rough and DOF) 
phases. 

We monitor also the density of AB pairs of parti- 
cles, pAs(t)- The A and B particles in each AB pair 
do not have to be immediately adjacent to each other. 
They can be separated by a stretch of empty space of 
arbitrary length. Similarly, pAA(t) and pss(t) are the 
densities of AA and BB pairs. It is easy to show that 
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p A (t) = PAA{t) + \pAB{t) and p B (t) = PBB{t) + \p AB {t). 
PAA{t) — Pb b (t) is preserved by the dynamics and equal 
to zero for nontilted surface configurations, paa distin- 
guishes between the DOF and the rough active phases. 
In the DOF stationary state, the system is active with 
perfect AB alternating order, i.e., p =fi and paa = 0. 
In the rough stationary state, the AB alternating order 
is broken, i.e., p ^ and paa ^ 0. 



III. THE DIRECTED ISING DYNAMIC 
PREROUGHENING TRANSITION 



In this section we present numerical evidence for the 
DI nature of the dynamic prcroughcning transition. Con- 
sider the f = line of the phase diagram (Fig. ||) where 
OBB type branching is excluded. Here, the configura- 
tions with perfect ABAB alternating DOF type order 
form a dynamically invariant subspace. 

In this subspace our model is almost identical to the 
BAW model with 2 offspring and dynamic branching [jl9| . 
There, each site may be empty or occupied by a particle 
of a single species. Those particles can hop to a nearest 
neighbor site with probability p or create two offspring on 
the nearest and next-nearest neighbor sites to the left or 
to the right with probability q = 1 — p (dynamic branch- 
ing). Two particles annihilate immediately when they 
happen to land on the same site. This BAW model ex- 
hibits an absorbing phase transition at p = 0.5105(7), 
which belongs to the DI universality class . 

There is an exact two-on-one mapping between the 
configurations in our model and those in the above BAW 
model. Simply label the particles in the latter as A and 
B (or B and A) alternatingly. The dynamic processes 
for the two models are virtually the same, except for one 
detail. In the BAW model, without A and B labels, par- 
ticles can always annihilate when they land on the same 
site. In our model only A and B pairs can annihilate, and 
AA and BB pairs repel each other. Hopping events are 
not affected by this, but OPB type branching processes 
are. Consider a configuration with an isolated AB near- 
est pair, like OOOABOOO where represents a vacant site. 
Suppose that the A particle is chosen to branch a pair 
of BA particles to the right. In our model, this attempt 
is rejected because it would result in two B particles on 
a single site. In the the BAW model this attempt is ac- 
cepted, and results in the annihilation of two particles. 
This difference between the two models does not affect 
the mod 2 conservation of the total particle number. So 
we expect a Dl-type absorbing phase transition along the 
f = line, but the OF phase must be more stable. 

The mapping between configurations of the BAW 
model and our model is lost outside the configurational 
subspace with perfect ABAB order. The DOF subspace 
is an attractor, however. It contains the global station- 
ary state for arbitrary initial configurations. Consider an 
arbitrary configuration, i.e., one with AA and BB pairs. 



For f = there is no mechanism to increase their to- 
tal numbers, Naa and Nbb- Hopping and OPB type 
branching decrease Naa and Nbb, by annihilation of 
AB pairs but never increase them. For example, con- 
sider a configuration like 0^40^4-8050. Hopping of the 
central A particle to the right induces the annihilation 
of an AB pair and the configuration becomes 0^40000-BO. 
Naa and Nbb decrease by one. The density of these 
pairs, paa — Pbb, decreases monotonically in time. In 
fact, they decay algebraically to zero along the entire 
f = axis (see section 5 and 6). The active steady state 
involves only configurations with perfect ABAB alter- 
nating order. 

To locate the critical point we perform defect dynamics 
type Monte Carlo simulations, in which one starts with a 
single nearest-neighbor pair of AB particles at the center 
of an empty lattice. Obviously, this initial configuration 
belongs to the DOF type invariant subspace. Time incre- 
ments by one unit after L single site updates (one Monte 
Carlo step) with L the lattice size. We measure the sur- 
vival probability P(t) (the probability that the system is 
still active at time t) and the number of kinks (particles) 
N(t) averaged over 5 x 10 5 independent runs up to 5000 
time steps. 

At criticality the long time limits of these two quanti- 
ties are governed by power laws with critical exponents 
S and 77 as P(t) ~ t~ s and N(t) ~ f f|. Precise esti- 
mates for the critical point and the critical exponents are 
obtained from a finite time analysis of 



■S(t) = 

v(t) = 



\n[P{t)/P{t/m)} 

lnm 
\n[N{t)/N{t/m)] 

lnm 



(1) 
(2) 



In Fig. [|, we plot these effective exponents against 1/t 
with m — 5 for several values of q. These plots bend up 
or down in time except at criticality. This leads to an 
estimate of the critical point, q c — 1.245(5) p2]], and also 
for the critical exponents, 



5 = 0.270(5) , ?7 = O.OI5(IO). 



(3) 



These values are slightly different from the standard DI 
values; 5 = 0.285 and 77 = 0.000. However, their sum is 
related to the steady-state properties via the generalized 
hyperscaling relation pi| and is in excellent agreement 
with the value of the DI universality class. 

Further evidence of the DI nature of the transition is 
obtained by monitoring the time evolution of the step 
density p{Sq, t, L) in stationary Monte Carlo simulations 
on a finite lattice of size L with Sq = q — q c . We start 
with a random configuration inside the invariant sub- 
space with initial kink densities pA (0) = /?b(0) = 1/4 and 
periodic boundary conditions. The step density p aver- 
aged over survived samples only should obey the scaling 
relation 



p(Sq, t, L) = b-M^pQbV^Sq, b- z t, b- l L). 



(4) 
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P is the order parameter exponent, u± the correlation 
length exponent in the spatial direction, z — v\\/v±_ the 
dynamic exponent p4|], u\\ the correlation length (relax- 
ation time) exponent in the time direction, and b an ar- 
bitrary scale factor. 

This scaling relation determines all scaling proper- 
ties of p. At the critical point, p decays algebraically; 
p ~ t~P/ u \\ for t <^tl with tx a characteristic time scale 
which diverges as L z . p scales at p ~ L^ l3 ^ u± for t > n,. 

Log-log plots of /o(0, L) versus t are shown in Fig. ||(a) 
at criticality for L = 2 5 ,...,2 10 . The step densities 
are averaged over 10 3 independent Monte Carlo runs at 
q = q c = 1.245. The inset shows the saturated values of 
the step density, denoted by po, as function of L. From 
these plots we obtain 

— = 0.282(5), —= 0.497(5), (5) 
"II v ± 

and z = 1.76(1) from their ratio. These results are in ex- 
cellent agreement with those of the DI universality class. 

The above analysis at the critical point gives only ra- 
tios of critical exponents. Their bare values can be ex- 
tracted from the off-critical behavior of p. The satu- 
rated value of the step density follows the scaling form 
p {Sq,L) = L-P/"±g(Sq L 1 /^), according to Eq. (g). 
The scaling function Q{x) becomes a constant in the 
x — > limit and scales as Q ~ x® in the x — > oo limit, 
because po ~ (Sq) 13 in the L — > oo limit. In Fig. 0(b), 
we show the log-log plots of poL 13 ^ 1 ^ against 5q L~^ ± . 
Using the values of q c obtained from the defect dynam- 
ics simulations and (3/v± from Eq. (g), we find that the 
data are best collapsed with v±_ = 1.78(5). Combining 
this value with those in Eq. (||), we obtain the critical 
exponents 

/3 = 0.88(3), v ± = 1.78(5), ^ = 3.12(5) . (6) 

These values agree well with the DI values and satisfy the 
generalized hyperscaling relation p3| , S + r] = (y± — 0)/v\\ 
very well. We conclude that our model has an absorbing 
type preroughening transition between the ordered flat 
and the DOF phase, which belongs to the DI universal- 
ity class. 

IV. CROSSOVER INTO THE ROUGH PHASE 

The PR phase transition is unstable with respect to 
OBB type branching, for all f > 0. We perform defect 
dynamics simulations at various values of r and q. The 
plots for the effective exponents, S(t) and r)(t), show only 
upward curvature. This indicates that the system is al- 
ways active for all f > 0. The size of the active region 
R(t) averaged over survived samples confirms this. It 
grows linearly in time so the spreading velocity of the ac- 
tive region is finite in the long time limit. In this active 
phase, the densities of AA and BB pairs, paa and pbb, 



are non-zero. This suggests that the steady-state surface 
is rough. We measure the surface width W(t, L) 

W{t 1 L) = l^ L j2{h l {t)~h{t)) 2 ^ (7) 

in Monte Carlo simulations on a finite lattice of size L. 
hi(t) is the height of the surface at site i, h(t) the average 
height, and (■ ■ •) the average over survived samples. In 
Fig. |[ we illustrate the typical behaviour of W(t, L) by 
the evolution of the surface at q = r = 1.0, using as initial 
condition • • • 0A0B0 ■ ■ ■ with system size L — 2 5 , . . . , 2 9 
and taking the average over 10 3 runs. The data satisfy 
the dynamic scaling form W{t,L) = L 2a f(t/L a ^) with 
the Edwards- Wilkinson roughness exponent a = 1/2 and 
the growth exponent f3 = 1/4. The active dynamic rough 
phase has the same scaling properties as the conventional 
equilibrium rough phase in one dimension. 

The crossover exponent y r at the DI critical point into 
the f direction must be relevant. y r is the scaling di- 
mension of the OBB type branching operator similar to 
Dt = X/v±_ which is the scaling dimension of the OPB 
type branching operator. y r is potentially an indepen- 
dent DI critical exponent (see Fig. |^). We obtain y r 
numerically by measuring the size of the active region 
R at various values of f along the q — q c line in defect 
dynamics type simulations. Consider the scaling rela- 
tion R(t,f) = bR(b~ ZDI t,b Vr r) with z DI = v\\/v±_ the DI 
dynamic exponent, and b an arbitrary scale factor. For 
b = t 1 ^ 01 , R takes the form R(t, f) = t x l Zm T{r W* DI ). 
The scaling function T(x) becomes a constant in the 
x — > limit, because R ~ t 1 /* 01 at the DI transi- 
tion point. In the x — > oo limit, T must scale as 
T ~ x^ Zol ~ 1 ^ Vr , because R grows linearly in time in 
the long time limit for f > 0. Therefore the asymp- 
totic value of the spreading velocity of the active region, 
v = limt^oo R/t, scales as v ~ f K with k = (zdi — X)/y r . 

In our simulations, we measure R up to times t = 
2 x 10 4 and average over 2 x 10 3 ~ 5 x 10 3 samples. 
The spreading velocity v is extracted in two ways. First, 
we fit R(t) to the form au + v^t in the time interval 
I0 k < t < 10 k+1 (k = 1,2,3,4). v k will converge to 
v in the asymptotic regime. Next, we define an effective 
velocity v eS (t) = (R(t + At) - R(t))/At with A< = 500, 
whose saturated value in the asymptotic regime is de- 
noted by Uoff- Fluctuations around the saturated value 
give an estimate of the statistical errors. In Fig. |[ we 
plot Vk and v e s for several values of f. The estimates for 
Vk with k > 3 merge into w c ff. This confirms that v e g is 
the asymptotic value of the spreading velocity v. From a 
power-law fit we estimate n — 0.64(1) and hence obtain 
the value of the crossover exponent, y r = 1.2(1). 

V. DENSITY DECAY DYNAMICS 

In section 3, we demonstrated the DI nature of the 
absorbing phase transition along r = 0. The active sta- 
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tionary state has perfect ABAB alternating order. The 
configurations with this perfect DOF type order form 
an invariant subspace. This subspace is an attractor, be- 
cause the density of AA and BB pairs, paa — Pbb, never 
increases in time in the absence of OBB type branching, 
along the f = line. In the conventional picture all 
asymptotic dynamic time scales depend only on the dy- 
namic exponent z and the stationary state exponents of 
the DI transition. Surprisingly, the latter is not true in 
our model. The step density 

P ~t- e (8) 

and the pair densities 

PAA = PBB ~ t-° AA , PAB ~ t- 0AB (9) 

decay in the large time limit with exponents that depend 
on the initial condition. Their values depend on whether 
the initial state has perfect AB alternating order or is 
rough, and whether this roughness is random or corre- 
lated. 

These densities decay as power laws everywhere along 
the f = line, with different exponents in the different 
phases. In this section we review first previously known 
results at point S of the phase diagram, see Fig. |[ and 
then generalize those results to the entire f — line. 

A. diffusion-limited pair annihilation: q — 

At point S (f = q = 0) all branching processes are 
disabled and the dynamics describes chemical reactions 
with diffusion-limited pair annihilation A + B — > 0. This 
process has been studied extensively p^-p8|. The parti- 
cles perform random walks, subject to an infinite on-site 
repulsion between the same species and AB pairs anni- 
hilate when they meet each other. A random initial con- 
figuration with equal initial densities, pa{0) — Pb(0) (a 
nontilted surface), decays to the absorbing OF "empty" 
state. 

Neglecting spatial correlations gives rise to the rate 
equation 

dpAjt) dp B (t) 

dt = dt =-pA{t)p B {t), (10) 

with p{t) ~ t^ 1 as solution. In the absence of on-site 
repulsion between the particles, it has been shown rig- 
orously that this mean-field behavior, p ~ t , holds 
in dimensions higher than 4 and that the particle den- 
sity decays as pit) ~ t~ e with 9 = d/4 in dimensions 
d < 4 P5fl . This contradicts the naive expectation that 
9 be equal to the inverse of the random walk dynamic 
exponent 9 = l/z IW = 1/2. 

The following scaling argument explains this result in- 
tuitively pq] . Let T>r be the difference in the number 
of A and B particles in a cube of size R d . For random 
(uncorrelated) initial configurations, Vr will be of order 



y/piO)R d . As the system evolves, each particle diffuses 
over a distance ~ i 1 / 2 ™ during time t (ignoring on-site 
repulsion). So after time tn ~ i? Zrw , all members of the 
minority species found a partner in the cube and have 
annihilated. This leaves the region occupied by the ma- 
jority species only. Therefore, the particle density decays 
as 

Pit) ~ 2V#W/«~ ~ VpW) t~ d/2z - ■ (11) 

This argument gives the correct value of the decay expo- 
nent 9 = d/2z rw = d/4 for d < 4. Several numerical sim- 
ulations |2^,|2^] found that the on-site repulsion between 
the same species does not alter the decay exponent 9. 

The same argument can be extended to the time evo- 
lution of the pair densities and the interparticle dis- 
tances [ p8| . We need them later in this section. The 
size of a domain occupied by one single species of parti- 
cles (e.g., a train of A's), Id, grows in time with the same 
exponent as the random walk radius, Id ~ i 1 / 2 ™. The 
AB pair density scales in d = 1 as pab ~ Id 1 because 
each domain with a single species of particles is bounded 
by two AB pairs. So pab decays with 9ab — l/^rw, he., 
decays faster than the particle density p; 9ab = 29. 

The AA and BB pair densities scale differently. They 
are equal for nontilted initial states. The pair densities 
add up to the total particle density, p = pab + Paa + Pbb 
(see section 2). Since pABit) decays faster than pit), 
Paa = Pbb must have the same asymptotic behavior as 
the particle density p, i.e., 9aa = 9 = l/2z IW in one 
dimension. 

The interparticle distances diverge in the long time 
limit. Define Iaa ^ab) as the average distance between 
the nearest neighboring particles of the same (different) 
species (see Fig. . These distances grow via pair annihi- 
lations of AB pairs. The A and B particle meet through 
random walk fluctuations at time intervals of order /Jg. 
Therefore, the particle density decays as dp(t)/dt ~ 

~PABit)/l Z AB, Which leads to l AB - i(2^w-l)/2z r 2 w _ ^3/8 

The sum of all interparticle distances adds up to the size 
of the system: 

ipAA + Pbb)Iaa + PabIab — 1- (12) 

The second term becomes negligible in the asymptotic 
limit, which yields Iaa ~ Paa ~ t x / 2z ™ . Iab diverges 
faster than Iaa but still slower than the interparticle dis- 
tance of the ordinary random walk problem. 

In summary, the dynamics at point S belongs to the 
random walk (diffusion) universality class. Starting with 
the nontilted random initial configurations, the density 
of AB pairs decays with the naive random walk dy- 
namic exponent, 9ab = l/^rw> but the density of AA 
pairs and the particle density decay much slower, with 
9 aa — 9 = l/2z rw . 

The above results assume that the initial condition is 
random. The exponents change when we modify the ini- 
tial state. The factor 2 appearing in the exponents 9 aa 



G 



and 9 originates from the random nature of the rough 
initial configurations discussed in the beginning of this 
subsection. Suppose it were correlated such that the ini- 
tial value for the surplus of particles of one species in a 
box of size R d scales with a different power, i.e., T)ji ~ 
[p(0)R d ] x . That leaves Qab = unchanged, but 

modifies the asymptotic behavior of the particle density 
and AA pair density to 9 = 9aa = (1 — x)/z l:w in one 
dimension. The exponents for the interparticle distances 
change accordingly. The nontilted random initial con- 
figurations correspond to x — 1/2 and the AB ordered 
initial configurations to x = 0. The tilted random initial 
configurations should correspond to x = 1. 

For initial configurations inside DOF subspace, paa 
and pbb are always equal to zero. Then the particle den- 
sity becomes equivalent to pab and decays with the "con- 
ventional" random walk dynamic exponent, p ~ t -1 / 2 ™. 

Let's now investigate how this dependence of the decay 
critical exponents on the initial states generalizes along 
the entire r = line, according to the same type of rea- 
soning. 



B. the absorbing phase: < q < q c 

At the absorbing side of the DI transition point, the 
asymptotic time scaling behavior remains the same as 
at point S except for one important detail. Clouds of 
particles take over the role of single particles. Consider 
an initial rough state with a low particle density. Each 
of these A and B particles broadens itself quickly into 
a small cloud of particles with a characteristic width £ 
via branching processes (see Fig. [?]). These clouds are 
created by OPB type branching, and therefore preserve 
ABAB alternating local order. This broadening is gov- 
erned by the DI type dynamics, like in defect dynamics 
with a single starting particle. The width of the clouds 
£ is finite and of the order of the DI correlation length. 
These clouds are well defined in the asymptotic limit be- 
cause the distances between them, the length scales, Iaa 
and Iab, diverge in time while £ remains finite. 

There are two topologically distinct types of clouds: 
the Ca clouds, nucleated from a single A particle and 
with j4's at both edges (AB ■ ■ ■ BA), and the Cb clouds 
(BA---AB). Clouds diffuse through hopping, branch- 
ing, and pair annihilation of bare particles. CaCb 
pairs can annihilate each other. Clouds do not branch. 
They could in principle, but a branching process like 
Ca —* GaCbCa, requires a collective sequence of mi- 
croscopic events involving many particles, and at length 
scales I > £ this does not happen. 

The clouds, at the DI length scale £, obey therefore the 
same dynamic rules as the diffusion-limited pair annihi- 
lation process of bare particles at point S, with renor- 
malized probabilities. The density of the clouds, p{t), 
and the cloud pair densities, paa (of CaCa cloud pairs) 
and pab (of CaCb pairs), must scale in the same way 



as particle densities at point S, i.e., p(t) ~ t 1 / 2z ™ : 
Paa (t) ~ t-V2»~ and pab (t) ~ i~ 1/2 ™. 

In numerical simulations we measure the bare particle 
densities. Those are related to the cloud densities as fol- 
lows. Each CaCa pair of clouds contains only one AA 
pair of bare particles, since each cloud consists out of a 
perfectly ordered ABAB train of particles. This implies 
that paa — Paa and 9aa — l/2z r w The number of AB 
pairs of bare particles in each cloud is proportional to its 
width £, so pab — £p(t) and 9ab = l/2£ r w The bare 
particle density is equal to the sum of all pair densities, 
P = Paa + Pbb + Pab- Therefore p scales with the slowest 
exponent, i.e., 6 = l/2z rw . 

The only difference with point S is that all three densi- 
ties decay with the same modified random walk exponent, 
~ t -1 / 22 ™. The random walk nature of the dynamics is 
completely obscured now. The random walk nature of 
the decay dynamics manifests itself only inside the sub- 
space with ABAB ordered initial configurations; there, 
9 = 1/zrw. 



C. the DOF active phase: q > q c 

In the active phase, solitons play the same role as the 
particle clouds do at the absorbing side of the DI phase 
transition. First, consider the limiting case q — oo (p — 
0,q = 1) where the hopping probability becomes neg- 
ligible with respect to OPB type branching. Any ran- 
dom initial configuration develops quickly by OPB type 
branching into fully-developed DOF domains separated 
by nearest-neighbor pairs of AA (step-up step-up) or 
BB (step-down step-down) particles as shown in Fig. 
These AA and BB pairs denoted by Sa and Sb in Fig. g 
are the topological excitations against the DOF phase 
and will be called the A- and B-type solitons. The den- 
sity of each soliton type is equal to the pair densities of 
bare particles, paa and pbb, respectively. 

At q — > oo, all activity is blocked inside each DOF 
domain, because any attempt of OPB type branching is 
rejected due to the infinite on-site repulsion between the 
same species and the fact that the DOF structure is close- 
packed. Only the solitons at the boundaries of the DOF 
domains are active. Their dynamics is basically identical 
to that of the bare particles at point S. Soliton diffusion 
is a second-order OPB process. Each soliton can hop to a 
next-nearest-neighbor site by applying OPB type branch- 
ing twice (Fig. ||). The solitons of the same species repel 
each other and SaSb soliton pairs can annihilate each 
other when they meet. There is no mechanism to create 
solitons. 

At finite q, the solitons broaden. Their width £ is 
the DI correlation length in the active phase. Just as 
the clouds in the absorbing phase, these broadened soli- 
tons must obey effective dynamics at length scales I > £ 
identical to those of the sharp solitons at q — > oo; i.e., 
identical to the particles at point S. The soliton density 
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decays with the same exponent 9 as that of the clouds at 
the other side of the transition. 

Since the soliton density is equal to the particle pair 
densities paa and pbb, we obtain paa ~ t~ 0AA with 
aa = l/2z rw - The particle density is finite in the steady 
active state. So p and pab do not decay as power laws, 
but remain nonzero. 



D. at the critical point: q — q c 

The decay dynamics at the DI critical point can be 
discussed equally well from the cloud or the soliton per- 
spective (the absorbing state or the active phase point of 
view). At the critical point the DI correlation length 
diverges in time as £ ~ f 1 / 2131 with the DI dynamic 
exponent zvi = v\\/v±_ — 1.76. This length diverges 
much faster than the typical distance between clouds 
(IcaCa ~ * 1/ ' 4 and IcaCb ~ t 3 / 8 ) (or solitons from the 
other point of view). This implies that the motions of the 
clouds (solitons) become correlated by DI critical fluctu- 
ations. Their diffusion is not driven by random walks 
with dynamic exponent z rw = 2, but by correlated ran- 
dom walks with the DI dynamic exponent zjyi. We expect 
this to be the only change in the scaling theory for the 
clouds (solitons). That means that we only need to re- 
place z rw by zni. The total density of clouds should scale 
as p(t) ~ t -1 / 2 *Di j the density of CaCa cloud pairs as 
PAAit) ~ t -1 ' 2 " 131 , and the density of CaCb cloud pairs 
as Pab (t) ~ t -1 /* 01 . The densities of the solitons (on the 
opposite side of the transition) scale identically. 

Next, we need to establish how these cloud and soli- 
ton densities are related to the bare particle densities. 
We did this for both, and the answer is the same, but 
the argument is easier from the clouds perspective and 
therefore we present only the former. 

The number of AA particle pairs is equal to the num- 
ber of CaCa soliton pairs (like above in subsection B). 
Therefore the AA particle and BB particle pair densities 
scale as PAAtt) = PBB(t) ~ f _1 / 2 *m. The total density 
of clouds scales as p(t) ~ ^-i/2«Dr j anc t therefore the av- 
erage distance between diverges as ln(t) ~ t 1 ' 2 * 01 . To 
find the total density of particles, we need to know how 
many clouds there are and how many particles each cloud 
carries. The width of the clouds, l w , diverges. For a sin- 
gle isolated cloud this happens with the DI correlation 
length, £(t) ~ t zm . This is much faster than the inter- 
cloud distance, Id- Therefore the width of the clouds is 
limited by the latter, l w = Id- The density of particles in- 
side each cloud scales as p p ~ lw P,vx (Eq. (|)), and their 
total number therefore as N p ~ ll~ 0/ " x ~ i(i-^/"x)/2a D i. 
The total density of steps scales as p ~ N p p. Putting all 
this together gives p ~ f^/P^Di) _ t -P/2u u g. 

The density of pab pairs is related to the other two 
by the relation p = IpAA + Pab and therefore must scale 
with the slowest power law, pab ~ t~ /3 / 2,y i . 



The density exponents are smaller by a factor of 2, 
compared to their asymptotic behavior starting from the 
AB ordered initial configurations (see section 3). This 
is the same factor of 2 found at point S and everywhere 
else along the f = line. This factor reflects the random 
roughness of the initial configurations. It changes for cor- 
related initial rough configurations in the same manner 
as discussed in subsection A for the S point, i.e., simply 
replace \ — > (1 — x). 

In conclusion, the above arguments predict that at 
the PR transition point the step density scales as p ~ 
j-a-aO/V^ii, that the AB step pair density scales with 
the same exponent, and that the AA and BB step den- 
sities scale as paa = pbb ~ t~^ x '' zm _ Zm ^ ^ anc j v ^ 
are directed Ising critical exponents, but x represents the 
correlations in the initial configuration. The initial state 
properties persist into the asymptotic scaling properties. 

VI. NUMERICAL SIMULATIONS FOR THE 
DENSITY DECAY 

The scaling theory in the previous section is heuristic, 
and certainly not exact. It is somewhat questionable in 
particular at the PR transition because we assume that 
the clouds (and solitons) remain valid concepts, while 
their widths actually like to diverge faster than allowed 
by the inter cloud (soliton) distances. 

To test these predictions, we perform Monte Carlo sim- 
ulations starting from a random initial state where parti- 
cles are distributed randomly on a lattice of size L = 2 15 
with initial densities pa = Pb = 1/4. We apply periodic 
boundary conditions. The time evolutions of the densi- 
ties p(t), pABit), and paa it) = Pbb (t) > are monitored up 
to t = 10 4 and averaged over 100 independent runs. A 
few simulations on a larger lattice of size L = 2 16 demon- 
strate that L = 2 15 is adequate to describe the scaling 
behavior up to time t = 10 4 . 

First, we test the density decay at point S (q = f = 0). 
In the absence of the on-site repulsion between the same 
species, the total particle density p should decay al- 
gebraically with exponent 9 = 1/4 pj|. It is impor- 
tant to confirm explicitly that the infinite on-site repul- 
sion between the same species in our model does not 
change this result. In Fig. |[ p is plotted against t on 
a log-log scale (the solid line). It seems that the den- 
sity decays slightly faster than £ -1 / 4 (the dashed line). 
Similar results were found previously |2q| . This is a 
correction-to-scaling effect. Insert the leading scalin g be - 
haviors of Iab ~ t 3 / 8 and pab ~ t- 1 ' 2 into Eq. @. 
This gives 2pAA x Iaa — (1 — at -1 / 8 ), and shows the 
presence of a generic 0(t -1 ' 8 ) type correction-to-scaling 
term. To isolate this term, we define an effective expo- 
nent 9{t) = -\n[p(t)/p(t/m)]/\nm with m = 8. The 
leading scaling exponent 9 is given by the limiting value 
of lirm^oo 9(t) and the correction-to-scaling behavior is 
contained in 9{t) - 9. The log-log plot of 0(f) - 1/4 
against 1/f is shown in the inset of Fig. ^. We find that 
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9(t) — 1/4 scales clearly as t 1 / 8 . This confirms that 
p(t) ~ <" 1 / 4 (1 + e»(^ 1 / 8 )). This correction-to-scaling 
term decays very slowly. Fitting p(t) to a simple power- 
law form therefore fails to produce the correct value of 
the leading exponent. We measured also the effective ex- 
ponents 9aa (t) and 9 ab (t) at point S. They approach 
1/4 and 1/2, respectively, with power-law corrections to 
scaling as well. 

The step density suffers from the same type of correc- 
tions to scaling everywhere along the f = line. Log-log 
plots of the density versus time are not straight lines. 
They curve a little. We analyze the data in the follow- 
ing manner. We construct estimates 9^ k ' (k — 1, 2, 3) for 
the exponent 9 by fitting the measured density p(t) to 
a power law in the time interval 10 k < t < 10 fc+1 (k — 
1,2,3). Approximants for 9^ A and 9^ B are constructed 
in the same way. As time increases the correction-to- 
scaling term contributes less and less, and the estimates 
should converge to the correct values of the leading decay 
exponents. 

In Fig. [l^, the estimates are presented and compared 
to the predicted values from the scaling theory in the pre- 
vious section. The step density exponent 9 is predicted 
to take the value of 9 = 1/4 below the transition and 
^{(3/v\\) ~ 0.141 at the critical point. Above the transi- 
tion, the density saturates to a finite value, i.e., 9 = 0. 
The estimates for 9 merge to the predicted value at the 
critical point q c = 1.245. Above the transition, 9^ be- 
come smaller as k increases, which is consistent with 
9 = 0. Below the transition, the convergence is slow 
(due to strong corrections to scaling) but compatible with 
61 = 1/4. 

The exponent 9 ab is predicted to take the value of 
9 ab = 1/2 at point S, 9ab = 1/4 along < q < q c , 
9ab = j{P/ v \\) — 0.141 at the transition point, and 
9 ab = in the DOF phase. The estimates, shown in 
Fig. [Io[ converge not as well as those for 9, but are still 
compatible with the theoretical predictions. The expo- 
nent 9 aa is expected to take the value of 9 aa = 1/2zdi — 
0.284 at the critical point and 9aa = 1/4 every where else. 
The data for this exponent converge slower than the other 
two, but they are compatible with the theoretical results 
as well. 



age with randomly placed kinks but in perfect up-down 
alternating order. 

The OF and DOF phases are both unstable with re- 
spect to the OBB type branching processes that break 
the up-down alternating order. There we find a rough 
stationary state with Edwards- Wilkinson type scaling be- 
haviour. The crossover exponent into the rough phase is 
determined numerically. 

The asymptotic long time behaviour of the kink density 
depends strongly on the statistical properties and correla- 
tions of the initial configurations. Information about the 
correlations in the rough initial state (controllable exper- 
imentally by sputtering for example) never gets lost. It 
obscures the random walk nature in the absorbing phase 
and the DI nature at criticality. We develop a scaling the- 
ory for the decay dynamics of the various kink densities 
and predict the values of the decay exponents. Numerical 
simulations confirm these. 

It is noteworthy to mention a different recent applica- 
tion of absorbing phase transitions to one-dimensional 
interface problems by Alon et al. fl3"of , which de- 
scribes the dynamic roughening phase transition from a 
smooth phase into a rough phase of Kardar-Parisi-Zhang 
type |31 1 . They considered a solid-on-solid type model 
where particles can adsorb at any site but desorption 
takes place only at existing steps. In other words, one 
can build mountains but is not allowed to dig new holes. 
This leads to the dynamic roughening phase transition at 
a finite value of the adsorption rate, which is triggered by 
the absorbing nature of the lowest level. This phase tran- 
sition belongs to the DP universality class. It may be in- 
teresting to introduce in our model a symmetry-breaking 
field between adsorption and desorption processes like in 
the above model. Generalization of our model in this 
direction is currently under study. 
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VII. CONCLUSIONS 



In this paper we presented a surface physics applica- 
tion of dynamic phase transitions in the directed Ising 
universality class. We introduced a model for the relax- 
ation of a one-dimensional interface starting from arbi- 
trary initial (rough) configurations, e.g., the straighten- 
ing of a step on a vicinal surface. This model undergoes 
a DI type dynamic preroughening transition between a 
perfectly ordered flat (OF) (absorbing) stationary state 
and a disordered flat (DOF) (active) stationary phase. 
The step becomes perfectly straight or straight in aver- 
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FIG. 1. Step hopping and branching processes of an up 
step (an A type particle). The letters s and d represent static 
and dynamic type branching, respectively. The dashed lines 
denote the surface before each event. Dynamic events near B 
particles are identical by mirror symmetry. 
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FIG. 2. Phase diagram of our model, f = r/p is the prob- 
ability for OBB type branching (creating local surface rough- 
ness) and q = q/p is the probability for OPB type branch- 
ing (creating disordered flat type local order), relative to the 
hopping probability p. y r and yr = l/v± are the scaling 
dimensions of the OBB and OPB operators, respectively. 
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FIG. 3. Semi-log plots of the effective exponents r\ and 5 
versus 1 /t for several values of q close to criticality. The data 
at our best estimate for the critical point q c are highlighted 
as thick lines. 
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FIG. 4. (a) Decay of the step density at q — q c . 
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tom. From the asymptotic slope of the curves we obtain 
/3/z/|| = 0.282(5). The broken line is of slope 0.282. The in- 
set shows the finite size scaling of the saturated step density. 
From the slope we obtain (3/v± = 0.497(5). The broken line 
has slope 0.497. (b) Scaling plot for log 10 [ j oo£ /3 ^ i ] against 
log 10 [(5gL 1 ^" i ] . Using the values of q~ c from defect dynamics 
simulations and (3/v± in (a), the best data collapse is obtained 
with u± = 1.78. 
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FIG. 5. The surface width in the rough phase at 
q = f = 1.0. The curves correspond to system sizes 
L — 2 5 , . . . , 2 9 . The width obeys the scaling form 
W = L 2a f(t/L a/fJ ) with the Edwards- Wilkinson exponents 
a = 1/2, /? = 1/4, and z — a/f3 = 2 as shown in the inset. 
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FIG. 6. Scaling of the spreading velocity v of the active 
region. Different symbols are used for vi (+), V2 (x), V3 (o), 
V4 (A), and v e g. (□). The Vk>3 merge into v c g., which con- 
firms that v e fi. is the asymptotic value for v. 



FIG. 9. The particle density p vs. t at point S. The dashed 
line has the expected slope —1/4. The inset shows that the 
effective exponent 6(t) approaches 1/4 with a power law cor- 
rection, i.e., 9(t) — 1/4 ~ i -1 / 8 . The dashed line in the inset 
has the predicted slope 1/8. 
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FIG. 7. Schematic typical surface configurations at the ab- 
sorbing side of the DI transition, < q < q c . 




FIG. 8. Schematic surface configurations at q = 1. They 
illustrate that hopping of the solitons requires a sequence of 
two branching processes. 
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FIG. 10. Estimates for the density exponents. The 
squares (□), up-triangles (A), and down-triangles (v) cor- 
respond to k = 1, 2, and 3, respectively. The predicted values 
of the exponents are represented by solid lines and circles. 
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